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Abstract. We examine here the constraints from the amount of relic density of neutralino 
dark matter and other experiments have on the SUSY parameter space for the mSUGRA model 
and for models with non-universal soft breaking at the GUT scale. In mSUGRA, the allowed 
amount of dark matter restricts the SUSY parameter space to a narrow band in toq — mi /2 
(except at very large tan/3). The Higgs mass and b ^ sj constraints produce a lower bound 
of TOi /2 ~ 300GeV and if the muon magnetic moment anomaly can be interpreted as a 3a 

deviation from the Standard Model, one also obtains an upper bound of mi /2 ~ 900GeV, 
making the SUSY spectrum well accessible to the LHC. The Bg —> fifi decay is seen to be 

accessible to the Tevatron Run2B with 15 fb~^ for tan/3 ~ 30. However, only parts of the 
spectrum will be accessible to the NLC if it’s energy is below SOOGeV. Non-universal soft 
breaking opens new regions of parameter space. Thus the mi /2 lower bound constraint of 
b ^ s'y and also the Higgs mass can be reduced greatly if the gluino mass is assumed larger 
at the GUT scale (allowing for a lighter gaugino spectrum), and non-universal Higgs soft 
breaking masses at the GUT scale can open new allowed regions at relatively low mi /2 and 
high tuq where dark matter detection cross sections may be increased by a factor of ten or 
more. 


1. Introduction 

Supersymmetry is a natural solution to the gauge hierarehy problem that sets in for the 
Standard Model (SM) at the TeV seale. Thus if one supersymmetrizes the Standard Model 
partiele speetra, one ean build a model going past the TeV scale yet consistent with all 
the successes of the Standard Model below 1 TeV. If one continues such models to yet 
higher energies, one find the remarkable grand unification of the three gauge coupling 
constants at the GUT scale Me = 2 x 10^® GeV, a result consistent with the LEP data 
at the percent level. Theoretical models which yield this unification arise naturally in 
supergravity (SUGRA) grand unification[[T|, and such models with R-parity invariance have 
the additional remarkable feature of predicting the existence of cold dark matter (CDM)[^, ^], 
the lightest neutralino x?, with a relic density amount comparable to what is observed. 

We discuss here some of the other consequences that might be expected of such SUGRA 
models. Already existing experiments have begun to restrict the SUSY parameter space 
significantly. Most significant of these are the amount of CDM, the Higgs mass bound, 
the 6 —S7 branching ratio, and (possibly) the muon anomaly. We start the discussion 
with the simplest model, mSUGRA, with universal soft breaking parameters, and discuss 
additional signatures of SUSY that might be seen at accelerators, i.e. the > /r+ + /i“ 
decay at the Tevatron, and the processes e+ + e“ — > + ff or e+ + e“ —> + X2 

future linear colliders (LC). We will then consider non-universal models (with non-universal 
gaugino masses and non universal Higgs soft breaking masses at Mg) to see how robust the 
mSUGRA predictions are, and where non-universal soft breaking signals might reside. 
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The mSUGRA models[j^, ^ depend on four extra parameters and one sign (and as such is 
the most predictive of the SUSY models). We take these to be mo (the universal scalar 
soft breaking mass at Mq), mi /2 (the universal gaugino masses as Mq), Aq (the cubic soft 
breaking mass at Mq), tan/3 =< H 2 > / < Hi > (at the electroweak scale), and the sign 
of the Higgs mixing parameter (which appears in the superpotential W as fcHiH 2 ). We 
examine the parameter range mo > 0, mi /2 <1 TeV (which corresponds to the LHC reach of 
rUg < 2.5 TeV), 2 < tan/3 < 55 and |Ao| < 4mi/2. 

mSUGRA makes predictions about two items involving dark matter: the neutralino- 
nucleus cross section being looked for by terrestrial detectors of dark matter from the Milky 
Way halo, and the mean amount of relic dark matter in the universe left over from the Big 
Bang. For detectors with heavy nuclei in their targets, the spin independent Xi-nucleus cross 
section dominates, which allows one to extract the x?—proton cross section o'^o_p. The basic 
quark diagrams here are the scattering through s-channel squarks, and t-channel CP even 
Higgs bosons, h and H. For the relic density analysis, one must calculate the neutralino 
annihilation cross section in the early universe which proceeds through s channel Z and 
Higgs poles {h, H, A) and t-channel sfermions. However, if a second particle becomes nearly 
degenerate with the one must include it in the annihilation channels which leads to the 
phenomena of co-annihilation. In SUGRA models, this accidental near degeneracy occurs 
naturally for the light stau, fi. One can see this qualitatively for low and intermediate tanbeta 
where one can solve the RGB analytically. Thus for the right selectron, cr, one finds 

mg^ = mg + 0.15mi/2 — sin‘^9wM^cos2l3, (1) 

M|o = 0 . 16 m 2/2 ( 2 ) 

For mo = 0, the two become degenerate at mi /2 = 350 GeV, and co-annihilation thus 
begins at mi /2 = 350 GeV (more precisely for the ri which is the lightest slepton). As 
mi /2 increases, mo must increase (to keep the fi heavier than the x?) and one gets narrow 
bands of allowed regions in the mo — mi /2 plane for each tan /3 and Aq. 

One starts the analysis at the Mq and uses the renormalization group equations (RGB) to 
obtain predictions at the electroweak scale. In carrying out these calculation, it is necessary 
to include a number of corrections and we list some of these here: (1) We use two loop gauge 
and one loop Yukawa RGB in running from Me to the electroweak scale Mew, and three loop 
QCD RGB below Mew for light quark contributions. (2) Two loop and pole mass corrections 
are included in the calculation of m/^. (3) One loop correction to irib and m.^ are included!^ 
(4) All stau-neutralino co-annihilation channels are included in the relic density calculation 
(Chargino-neutralino co-annihilation does not occur for mi /2 < 1 TeV.) Barge 
tan /3 NBO SUSY corrections to b ^ s'y are included[^, [^, PT] ]. We do not include Yukawa 
unifications or proton decay constraints as these depend sensitively on post-GUT physics, 
about which little is known. 

We use the following experimental input to constrain the SUSY parameter space: (1) 
Global fits for the CMB and other astronomical measurements now restrict the amount of 
CDM considerably [|T^ . We use here the range: 

0.07 < klcDuh'^ < 0.21 (3) 

The MAP data, due out perhaps early next year should be able to restrict this range 
considerably. 

(2) The BBP lower bound ruh >114 GeV[jT^, is a significant constraint for tan /3 ~ 30, 
and in fact an increase of three to five GeV would cover most of the parameter space. 
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Figure 1. Example of a leading contribution to the decay Bs fifx. Each vertex with a dot 
has a factor of tanbeta so that the diagram is proportional to tan (3^. 
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Figure 2. Expected sensitivity to Bg p/i of the CDF detector as a function of luminosity. 
The solid curve is a conservative estimate, and the dotted curve gives a limiting sensitivity. [20] 


Unfortunately, however, the theoretieal evaluation of ruh still has an error of ~ (2 — 3)GeV, 
and so we interpret the LEP bound to mean > 111 GeV. 

(3) The CLEG data for the decay b —> 57 [p^ has both systematic and theoretical error, 
and so we use a relatively broad range for the branching ratio around the CEEO central data: 

1.8 X 10“^ < B{B Xs'j) < 4.5 X 10“^ (4) 


The b ^ S'j rate produces a significant constraint for large tanbeta. 

(4) Shortly after this conference, the Brookhaven E821 experiment published new results 
of their measurement of the muon anomaly [[T^ which reduced the statistical error by a factor 
of 2 (and the systematic error somewhat). Also new more accurate data from CMD-2, BES, 
AEEPH and CEEO have allowed a more accurate determination of the SM contribution to a^, 
and there have been two new evaluations of this [|T^ [1^ (See also talk by T. Teubner at this 
conference). If the e~^ — e~ data is used to calculate the SM contribution, both groups find a 
3cr deviation between experiment and the SM prediction e.g. [0 

Aa^ = 33.9(11.2) x 10"^° 


( 5 ) 































Dark matter and Consequences of SUSY 


4 



Figure 3. mSUGRA predicted values of B(Bs —s- /i/i) for different values of tanbeta for 
Aq = 0 and mo = 300. 


On the other hand if the tan data of ALEPH and CLEO are used (with appropriate CVC 
breakdown eorreetions) the deviation from the SM is reduced to only 1.6cr and the two 
evaluations are statistically inconsistent The matter remains unclear as to which result is 
correct. In order to see the significance of Eq. (5), we will here use this result, and assume a 
2cr lower bound on delta exists and that this is attributable to SUSY. SUSY, in fact makes 

a significant contribution to a^, and if were ~ 10“^°, the squark and gluino mass 

spectrum would be pushed into the TeV domain. 

The combination of the m/i(for low tan/3 ) and 6 —sy (for high tan/3) constraints 
produces a lower bound on mi /2 over the entire parameter space of mi /2 ~(300-400)GeV, 
and consequently m^o ~(120 - 160)GeV. This means that most of the parmaeter space is 
in the fi — x? co-annihilation domain, and hence in order to satisfy the CDM amount, uiq 
is approximately determined by mi/2 (for fixed tanbeta, Aq). The sign of /i and Aa^ are 
correlated [p^ P^ , and so the assumption of Eq. (5) means that /r > 0. 


3. SUSY Signals At The Tevatron and Linear Colliders 


Before examining in more detail the effects of existing experiments on the SUSY parameter 
space, we discuss the Bg /x/i decay at the Tevatron and also examine SUSY signals that 
might survive at a linear collider. 

The Bg icfj. opens an additional window at the Tevatron for investigating the 
mSUGRA parameter space[|2^, ^ The SM predicts a very small branching ratio 
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Figure 4. mSUGRA allowed region for tan /3 = 10, Aq = 0, fj, > 0. The shaded upper right 
region is forbidden by the bound at the 2cr level if Eq. (5) is valid. The dot-dash lines are 
for the LC Xi ~ X 2 signal for the 500GeV machine (left line) and SOOGeV machine (right 
line). (The LC is sensitive to regions to left of the line.) The curved solid lines are for the LC 
fi — fi signal (the lower one for the 500 GeV machine and the higher one for the SOOGeV 
machine). The short vertical lines are for DM detection cross sections, cr^o_p = 5 x 10“®pb 
(left line) and 1 x 10“®pb (right). 


of B[i?s ^ /i/i] = (3.1 ± 1.4) X 10“®. Further, the SUSY eontribution ean beeome quite 
large for large tanbeta, as the leading diagrams (an example is given in Fig. 1) grow like 
tan/?^ and henee the branehing ratio grows as tan/?®. A set of outs have been obtained 
eliminating baekground (mostly gluon splitting g ^ bb and fakes) [ p| ] leading to a sensitivity 
of B[i ?5 ^ fifi] ~ 6 X 10“® for 15fb~^ /deteetor. Fig. 2 shows the sensitivity (for a single 
detector) as a function of luminosity, and Fig. 3 shows the expected mSUGRA branching 
ratio for different tan (3 for the example of Aq = 0, mo = SOOGeV. One see that the Tevatron 
will be sensitive to this decay for tan (3 ~ 30. 

It is interesting to examine what possible SUSY signals of mSUGRA a linear collider 
might see. We consider here two possibilities: ^/s = 500GeV and a/s = SOOGeV. We’ve seen 
above that the ruh and b ^ s'y constraint already mean that mi /2 > (350 — 400) GeV, and 
so for mSUGRA this means that gluinos and squarks would generally be beyond the reach of 
such machines (as well as selectrons and smuons for a large part of the parameter space). The 
most favorable SUSY signals are then 

e+ + e- ^ X 2 + X? ^ (/+ + 1 - + X?) + X? (6) 

e+ + e" ^ f/ + rf ^ (r + x?) + (r + x?) 

where or l~ means any charged lepton. Since for mSUGRA one has that m^o ~ 2m^o, the 
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Figure 5. The same as Fig. 4 for tan [3 = 40, Aq = 0, the dashed lines giving contours of 
Bs uu branching ratios. The lower short vertical line is cr^o_p = 3 x 10“®pb and the upper 
one is 1 X 10“®pb 


mass reach for these particle are 

l/2m^o ~ ~ 165(265) (7) 

TUf^ ~ 250(400) GeV for = 500(800) (8) 

There are a number of SM backgrounds that have to be considered to see if a clean 
signal remains. For example, WW production with decay into tau final states can be 
controlled by polarizing the beams. Cuts to eliminate other possible backgrounds are under 
consideration . 

4. The mSUGRA Parameter Space 

We now summarize the effects of the constraints from irih, b sj, dark matter density and 
on the mSUGRA parameter space, and also what might be expected from the Bg —> /x/i 
observation at the Tevatron and linear collider (LC) signals. Fig. 4 exhibits the parameter 
space for tan (3 = 10, Aq = 0. One sees that the bound (if valid) combined with the ruh 
bound and the relic density constraint leaves very little parameter space at low tan f3. Either of 
the two linear collider signals could scan the full remaining space. The dark matter detection 
cross sections are of size that would be observable by the planned future DM detectors. Figs. 
5,6,7 for tan/? = 40, Aq = 0, 2mi/2 and -2mi/2 exhibit a larger allowed parameter space. 
The Tevatron’s Bg /i/i covers the full allowed region for Aq = 0 and - 2mi/2, and about 
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Figure 6. Same as Fig. 5 for tan [3 = 40, Aq = 2mii2, /r > 0. 


half the space for ^40 = 2mi/2. In each case, the NLC at 500GeV can only cover a part of 
the parameter space. We note that the stau-stau LC signal extends to high mi /2 and relatively 
low mo, while the X? - X2 signal can cover large ttiq but limited mi/ 2 . In view of the nature 
of the allowed DM channel, the former is of more significance. At very high tan (5, a bulge in 
the DM allowed channel develops at low mi /2 due to the fact the heavy Higgs {A,H) become 
light, allowing a rapid early universe annihilation through the A and H s-channel diagrams. 
This is shown most dramatically in Fig. 8 for tan/3 = 55, Aq = 0. Here even the 800 GeV 
LC cannot cover the full parameter space, though the Tevatron signal of Bg would be 

observable over the full parameter space. 

5. Non-Universal Models 

We examine next some SUGRA models with non-universal soft breaking to see what aspects 
of the results of Sec. 4 are maintained. We consider specifically the case of non-universal 
gaugino masses and non-universal Higgs soft breaking masses at MG. Two striking effects 
in mSUGRA are the Xi ~ A co-annihilation channel which leads to a narrow band in the 
nr 1/2 — mo plane of allowed relic density, and the fact that the combined effects of 6 —> sy 

and the mh bounds require m^o ~ 120 (i.e. mi /2 ~ SOOGeV). 

The possibility of non-universal gaugino masses at Mg can relax significantly the 
constraints of b ^ s'f and mh- Thus if the gluino mass is increased, it effects the stop 
mass which has major effects on both 5 —> sy and mh- For example, a gluino mass which is 
twice the universal value at Mg reduces the lower bound on m^o to 75GeV (mi /2 ~ lOOGeV) 
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Figure 7. Same as Fig. 5 for tan j3 = 40, Aq = —2mi/2, /i > 0 


and also m-± to ~ ISOGeV, making these partieles more aeeessible. On the other hand, the 

near degeneraey between fi and x? remains, sinee this depends mainly on the U(l) gaugino 
mass rfii- Non-universality in the Higgs masses at Mq, i. e. = ml{l -|- 51 ^ 2)7 can 

also produee new effeets. While this introduees two additional parameters into the model, 
one ean understand their effeets sinee eontrols mueh of the physies. Thus if /i^ deereases, 
the Higgsino eontent of the neutralino inereases and this has two effeets: it inereases the 
x\-x\-Z eoupling and also To see what effeets oeeur qualitatively, we note that for 

small and intermediate tan /3 the RGB ean be solved analytieally and give 

= (/^^)mSUGRA + — 1)[—1/2(1 + Dq)52 + (9) 

where t = tan/3, Dq = 1 — (mt/200sm/3)^ = 0.25. A positive 82 then deereases jjf while 
jjf is relatively insensitive to (5i. Fig. 9 exhibits the allowed relie density regions for the ease 
of 82 = 1, tan (3 = 40, Aq = mi/ 2 . One sees that a new region of allowed relie density arises 
for small mi /2 and large mo from the inereased annihilation through the Z s-ehannel. (The 
usual fi — Xi eo-annihilation region is present also.) This region is mostly unobservable at 
a 500GeV NBC. Fig. 10 shows the eorresponding DM eross seetions. The Z-ehannel region 
ean inerease cr^o_p by a faetor of 10 or more. 

6. Conclusions and Summary 

We have diseussed here the eurrent restrietions on the SUSY parameter spaee from existing 
experiments and what may be obtained from future experiments at the Tevatron (from the 
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Figure 8. Same as Fig. 5 for tan/3 = 55, Aq = 0, /x > 0 and mt = 175GeV, mt = 4.25. 
(Note that the results are sensitive to the exact values of mt and mt, used 



Bs /x/i decay) and from linear colliders (using fi — fi and Xi — X 2 signals). We have 
examined mSUGRA models and also non-universal SUGRA models. 

For mSUGRA one finds the fi — x? co-annihilation limits the allowed parameter space 
to a narrow band except for large tan [3. The combined b ^ sj and rrih bounds then requires 
m^o ~120GeV, forbidding a light neutralino. The Tevatron with 15 fb~^/detector luminosity 

could scan almost the entire parameter space for tan {3 40 (if ~ 10 X 10^°) using 

the Bg —/x/i decay. A linear collider at SOOGeV can scan almost all the parameter space, 
though a 500 GeV machine would miss much of the parameter space for large tan f3. Dark 
matter detection cross sections are all within the range of sensitivity for planned detectors (i.e. 
~ 10-^0 pb). 

Non-universal models introduce new phenomena. Thus a non-universal gluino mass can 
greatly weaken the lower bound on m^o due to 6 —57 and m/j. An increased non-universal 
H 2 Higgs mass at Mq can give rise to a new region of allowed relic density for low mi /2 
and large mo from rapid Z- channel annihilation, while maintaining the usual fi — 7 ° co¬ 
annihilation channel. The dark matter cross sections can increase by a factor of 10 or more in 
the new Z-channel region. 
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Figure 9. Same as Fig. 5 for <52 = 1, tan [3 = 40, Aq = miji, q > 0- The shaded region for 
large toq low mi /2 is due to the increased annihilation through the Z-channel. The lower 
stau-neutralino narrow band is essentially unchanged. 



Figure 10. Dark matter detection cross section for tan/3 = 40, = rni/ 2 , q > 0, 62 = 1- 

The upper dashed band is the new region due to the increased Z-channel annihilation. The 
lower line us the upper value from the narrow fi — co-annihilation channel. The short 
vertical lines are the 2 a bound from the anomaly. 
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